In [114]:
# Standard setup block for running Python code
from __future__ import division, print_function
import os
if os.path.split(os.getcwd())[-1] == "Lab notebooks":
os.chdir("../../")
print("Moved to experiment root directory")
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
plt.style.use("Config/plotstyle.mplstyle")
from Modules import processing as pr, plotting as pl
%matplotlib inline
Got to the lab about 1:30 PM.
Homed turbine axis.
Wiped corrosion from turbine mounting frame. No visible corrosion on turbine parts.
Tank depth is about 2.45 m.
Homed tow axis.
Checked outer guy wire tension. Looks good.
Homed y- and z- axes.
z-axis missed some steps on upward motion. Will grease ball screw and nut then re-home.
3:22 PM -- Greased up z-axis ball screw and nut. Binding is gone.
3:33 PM -- Going to do a shakedown run (no. 20) at 0.8 m/s $\lambda=3.4$ to see how we're doing for finding $\lambda_0$. Obtained $C_P = 0.36$
3:49 PM -- Shakedown run 21, $U_\infty=0.9$ m/s, $\lambda=3.3$, $C_P=0.37$
3:53 PM -- Shakedown run 22, $U_\infty=1.0$ m/s, $\lambda=3.3$, $C_P=0.37$ -- looks like we may have found $Re_0$.
3:57 PM -- Shakedown run 23, $U_\infty=0.3$ m/s, $\lambda=3.3$, $C_P=0.23$
4:06 PM -- Shakedown run 24, $U_\infty=0.7$ m/s, $\lambda=5.6$, $C_P=-0.04$ -- carriage oscillates a bit for the first 5 m or so. std_tow_speed = 0.015 m/s.
4:20 PM -- Shakedown run 25, $U_\infty = 0.8$, $\lambda = 3.25$, $C_P = 0.37$ -- could make tsr range from 0.25--5.5 with increments of 0.25.
4:48 PM -- Added a cup of seeding particles to tank. Did shakedown run at 0.8 m/s, 3.25 tsr. Vectrino SNR went up to about 15 dB.
4:58 PM -- Did shakedown run 27. $\lambda=5.5$, $U_\infty=0.8$. Carriage oscillates a bit at mid tank. $C_P=0.04$.
5:02 PM -- Did shakedown run 28. $\lambda=5.0$, $U_\infty=0.8$, $C_P=0.15$. std_tow_speed = 0.026. Significant carriage oscillation due to belt resonance.
7:44 PM -- Turbine zero angle bracket got caught on curved jaw coupling and got bent pretty badly. Will try to turn turbine on/off each day at zero angle to keep things repeatable.
9:37 PM -- Calling it quits for tonight. Pretty close to knowing TSR values for performance curve, but want to play around with some different options. Settling runs will be done tomorrow.
In [115]:
runs = [27, 28] + range(32, 74)
#runs = [34, 44, 48, 49, 52] + range(57, 66) + [67, 68, 69]
cp = []
tsr = []
for n in runs:
r = pr.Run("shakedown", n)
cp.append(r.mean_cp)
tsr.append(r.mean_tsr)
df = pd.DataFrame()
df["run"] = runs
df["tsr"] = tsr
df["cp"] = cp
print(df)
plt.figure(figsize=(10,6))
plt.plot(tsr, cp, "o")
plt.xlim((0, 5.75))
plt.ylim((0, 0.4))
plt.xlabel("$\lambda$")
plt.ylabel("$C_P$")
plt.show()
In [68]:
3.3*pr.R*.06*0.8/1e-6
Out[68]:
In [90]:
pr.process_latest_run("shakedown")
In [35]:
pr.process_run("shakedown", 37)
Out[35]:
In [ ]: